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New issues for Numerical Stochastic Perturbation Theory 
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First attempts in the application of Numerical Stochastic Perturbation Theory (NSPT) to the problem of 
pushing one loop further the computation of SU(3) (SU(2)) pertubative /3 function (in different schemes) are 
reviewed and the relevance of such a computation is discussed. Other issues include the proposal of a different 
strategy for gauge-fixed NSPT computations in lattice QCD. 
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1. Perturbative /3-functions by NSPT. 

Perturbative calculations in lattice field theo- 
ries are necessary to connect results in different 
renormalisation schemes. In asymptotically free 
theories particulary useful are matchings to those 
schemes in which perturbative computations at 
high energies are more easily achieved (first of all 
MS). It is well known that the bare definition of 
the coupling in lattice QCD is often not the best 
one: other definitions are useful, for instance, to 
relate the perturbative and the non perturbative 
regions in a smoother way. 

In the case of QCD with no (or massless) 
quarks, or in any other pure gauge theory, just 
one free parameter is present and every coupling 
constant (in any scheme) can be expanded as a 
power series in any other. The first coefficient 
of such an expansion determines the relation be- 
tween the physical scales (A parameter), the oth- 
ers are related to combinations of the respective 
(non universal) coefhcients of the /3-function. 

Aim of the project is to obtain higher order 
expansions of convenient couplings in the lattice 
scheme, thus obtaining informations on higher or- 
der /3-function coefficients. The natural goal is 
/Ja in the lattice scheme, which would result for 
example in a more precise correction to asymp- 
totic scaling and in pushing one loop further the 
matching to MS. 



2. NSPT for lattice QCD 

Unfortunately perturbation theory presents 
many difficulties in the lattice regularisation. To 
bypass such difficulties NSPT was introduced, 
which proved successful in performing high loops 
gauge invariant computations in lattice QCD [Q. 

The idea which makes NSPT possible is the fol- 
lowing: The Langevin equation provides a quan- 
tisation scheme for field theories, and also a tool 
for numerical computations of the Functional In- 
tegral. It may be introduced by giving the fields 
a dependence on a new parameter (stochastic 
time), [/(/ijx) — > Urj{^,,x,t), x £ ivit) is an 
appropriate random process), and then imposing 
the fields to obey the Langevin equation, 

d_ 

dt 



f/, = (-*V5[[/,]|,,,^,-zr;)C/, 
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which determines the evolution of the fields in 
stochastic time. One then notices in the Langevin 
equation the explicit dependence on g through the 
action S; as a result the solution f/^ depends on g 
as well and can be formally expanded as a power 
series in g. This produces a cascade of stochastic 
equations, of which only the first depends explic- 
itily on the noise rj and each one depends only on 
the lower order fields. 

3. Symanzik's analysis 

In order to get a continuum-limit result from 
perturbative calculations on the lattice it is nee- 
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essary to recognize in the data the dependence 
on the cutoff. Typically one computes a coupling 
which is defined at a length scale L on a lattice 
of size / = L/a, thus obtaining a generic pertur- 
bative coefficient rn{I); then Symanzik's analysis 
suggests an asymptotic expansion for m{I) in the 

form: m(/) = J2'^=o Y.k=Q ^^^"^T^r^, where the 
maximum power of the log is fixed by the order of 
the computation. It is actually this analysis that 
sets the required precision of the computation. 

4. NSPT for the Schroedinger Functional 
scheme 

One particulary interesting renormalisation 
scheme is introduced via the running cou- 
pling g'^{L) introduced in |^ and based on the 
Schroedinger Functional with an induced back- 
ground field, obtained by fixing boundary condi- 
tions in the t = and t — L time-slices. The cou- 
pling g^{L) is defined through the rensponse to a 
change in the boundary conditions Through 
a finite size scaling technique |^ the ALPHA Col- 
laboration has been able to compute this coupling 
on a huge momentum range up to the perturba- 
tive domain in which contact can be made with 
other schemes such as MS. Beeing this scheme a 
successful one with many respects and having per- 
turbative computations of this coupling already 
been performed at two loop level (so that cross- 
checks are possible), it was a natural candidate 
for our analysis. 

It is however the first time that NSPT is ap- 
plied starting from a vacuum which is not the 
trivial one. In fact the fields have to be rep- 
resented as perturbative expansion around the 
classical solution induced by the background field 
= ^+J2k>o U^''^ This poses no problem 
in principle, but it will be interesting to see which 
kind of effects such a modification brings to the 
precision of the results. As far as the observable 
is concerned, it is known to exhibit quite strong 
fluctuations even at non-perturbative level. A 
feature that immediately emerges is that it is ef- 
fectively defined only on the two boundary time 
slices, and not on the bulk. 

As it has been soon recognized NSPT in gauge 
field theories is hopeless without introducing any 



kind of stochastic gauge fixing (see below). In 
fact, even if this is not necessary in principle, di- 
verging fluctuation will soon dominate without it 
[^,|| . A natural way of doing it is to apply (after 
every Langevin evolution step) a gauge trasfor- 
mation UJ^^in) = e^t") f/^(n) e-'"("+'^) where 

= -f E^^^(t^M(") - t^^("))traceless ^nd 
(0 < a < 1). It is easy to verify that such a 
trasformation attracts the system towards config- 
urations where the norms of the fields are under 
control. One can also apply a slightly modified 
version of the Landau gauge prescription (because 
of the presence of a non trivial vacuum) , however 
this produces no significant difference in the re- 
sults. 

Unfortunately the results for the observable 
defining g^{L) are extremely noisy, preventing 
any Symanzyk analysis to get a continuum limit. 
Perturbative calculations had been carried on an- 
alytically in 1^ with very high precision to the 
second non trivial order. These known results 
are reproduced, for both SU{2) and SU{3), and 
it would be easy to go to higher orders. However 
the errorbars makes these results useless in the 
context they are interesting. 

Two considerations are in order: first the 
imprecision is due essentially to the observable 
and not to the new vacuum configuration upon 
which perturbations are made, nor to the differ- 
ent choice of gauge fixing one is forced to take. 
In fact simulations of the plaquette have the same 
precision we were used to. This suggests that this 
approach could be more profitably applied to ob- 
servables build up with Wilson loop, leaving still 
the freedom of choice of other conditions. Sec- 
ondly we convinced ourself that a great help could 
come from a more efficient procedure of gauge fix- 
ing 0. This is the subject of what follow. 

5. Fadeev & Popov in NSPT 

We now propose a new formalism which is even 
farer from the original spirit ("Perturbation 
Theory without Gauge Fixing"), but implements 
the Fadeev & Popov mechanism in a (perturba- 
tive) Monte Carlo. Such a formalism consists in 
the application to NSPT of a technique which is 
well known in the standard Langevin equation 
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and is in principle a solution for NSPT "dynam- 
ical fermions" as well. 

Consider any problem of the class det M 
(fermions, F& P). In principle e~"^detM ~ 
g_(s_Tr log M) ^ solution in the standard 
Langcvin equation has been known for a long 
timeiconsider U{t + 1) = e~^'^U{t) where fi = 
€[diSg - miMki^d^Mir,^] and ^ {k multi in- 
dex) are Gaussian as well with = 1. Again 
to 0(e) what matters in Fokker-Plank equation 
is (average on ^) {fi) = edi[Sg - Tr log M] + ^/erj 
(where 77 arc Gaussian as well with (77^) = 2) 
which means for the probability density of the 
fields P ~ e-('5-T'-iogM) ^ g-'^detM 

Instead of the fcrmion determinant consider the 
functional integral which is relevant for PT in 
covariant gauges (i.e. Fadeev & Popov mech- 
anism) (O) = 1/Z J DUe-'^^^+^<^f^AFp[U]0[U] 
where Sg = gauge Wilson action, Sgf is the gauge 
fixing action, and Ai?p[?7] the Fadeev Popov de- 
terminant. Sgf is in fact the action in charge of 
"killing the divergences" in Langevin equation. 

It is necessary to compute the contributions to 
the equation of motion coming both from Sgf and 
from /S.pp The first gives no problem. As far 
as the second is concerned notice that of course 
nobody implements A^? p via A pp = det L = 
gTr log The resulting action would be strongly 
non local. Instead deti is in standard applica- 
tions obtained by a functional integral over grass- 
mann fields (ghosts). Shadow of the non lo- 
cality of the would be Sfp[U] is in /j = 
e[di{Sg + Sgf)-^[C'<L-^diL^]] but the crucial ob- 
servation is now that L'^ = Lq^ — gL^^LiLQ^ + 
g^iLo'LiLo'L^L^' - L^'L^L^') + Oig^). The 
message is that has a simple recursive form 
easy to guess also at this order; only Lq ^ is a 
proper inverse, which is well known and easy in 
Fourier space and it does not depend on the fields. 

6. Conclusions 

NSPT proved successful in performing high 
loop computation in lattice QCD. The preci- 
sion demostrated in the computation of different 
quantities show that the task of pushing one loop 
further the evaluation of SU{3) {SU{2)) jS func- 
tion is attainable. 



However the choice of the scheme based on the 
Schroedinger Functional (which has many desir- 
able features) proved unlucky for our approach. 
The direction to pursue is now double. On one 
side one can examine different schemes of renor- 
malization (for instance based on Wilson loops), 
where NSPT calculations are more precise, and 
determine in this way the next coefficient of the 
(3 function. On the other side one can look at the 
possibility of implementing the Fadeev & Popov 
gauge fixing in NSPT. One notes that contribu- 
tions to L"'^ arc out of a recursive structure and 
depend only on one inverse L^^, which is easy 
to manage in Fourier space. L^^ does not de- 
pend on the fields. Because of the former observa- 
tion, the actual implementation of the proposed 
method depends strongly on the architecture of 
the computer: without a good FFT it is hopeless. 
Same sort of remarks apply to including dynam- 
ical fermions contributions as well. 
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